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مقاله علمي - پژوهشی 


بررسی دقت برخی از مدل‌ها در برآورد ضرایب آزمون تولید گاز در سیلاژ ذرت 


خلیل زابلی"۱- سارا کلوندی"- مصطفی ملکی! 


تاریخ دریافت: ۱۳۹۸/۰۲/۲۳ 


چکیده 


ه‌منظور پیش‌بینی کینتیک تخمیر شکمبه‌ای سیلام ذرت در زمان‌های ضفره ۳۰ و ۶۰ روز پس از سیلو کردن, از داد چهار مدل غیر خی امنفاه 
شد. برای این منظور. داده‌های تولید گاز در طول ۱۴۴ ساعت انکوباسیون با استفاده از مدل‌های نمایی (1726۳) لجستیک (100)» گومپرتز (/007) و 
فیزو (:12) برازش شدند. نکویی برازش مدل‌ها با استفاده از آماره‌های ميانگین مربعات خطا (4(۷8815 ضریب تعیین (4182 انحراف مطلق میانگین 
باقیمانده (18۷/۸) و میانگین درصد خطا (1۷1۳17) انجام شد. از آزمون‌های دوربین -واتسون (121۷): آزمون اجرا و رگرسیون خطی به‌منظور بررسی 
دقت مدل‌ها استفاده شد. نتایج نشان داد مدل‌ها از نظر پیش‌بینی پتانسیل تولید گاز (۸) تفاوت معنی‌داری با هم نداشتند. مدل "126 دارای بیشترین 
مقدار ۱455 (۲۹/۵۰) و کمترین مقدار 182 (۰/۹۶۱) در بین مدل‌ها بود. مقدار 131۸0 در مدل 0034 و ۳2۳ کمترین (به‌ترتیب ۲/۵۹۱ و ۲/۸۷۹) 
و در مدل 726۳ بیشترین (۳/۸۰۷) مقدار بود (0>۰/۰۵). در مدل ۳26۳ مقدار 1۷۳ ۵/۵۲۷ به‌دست آمد که تفاوت معنی‌داری با سایر مدل‌ها داشت 
(۰/۰۵>). در مدل 326۳ آماره 12۷۷ به عدد صفر نزدیک‌تر بود (۰/۳۹۲) که نشان‌دهنده ضعف مدل ۳26۳ در برازش داده‌ها در مقایسه با سایر مدل‌ها 
بود. بررسی رابطه رگرسیون خطی نشان داد که مدل‌های 211 و 1606 پیش‌بینی بهتری از پروفیل تولید گاز داشتند. به‌طور کلی نتایج نشان داد مدل - 
های 001۷ و ۳21 جهت پیش‌بینی کینتیک تخمیر شکمبه‌ای سیلاژ ذرت از دقت بیشتری در مقایسه با مدل 1726۳ برخوردار بودند. 


وازه‌های کلیدی: روش آزمون تولید گاز» کینتیک تخمیر شکمبه» مدل‌های غیر خطی و نکویی برازش 


مقدمه 

روش‌های متعددی برای ارزشیابی مواد خوراکی در تغذیه 
تفتکوار کتند کان. مود فارند کد.فر کنام از آنها دازای مزابا و مسایی 
هستند. استفاده از روش درون‌تنی" برای تعیین قابلیت هضم» دارای 
ارزش زیادی است. اما این روش پر هزینه و زمان‌بر بوده و نیاز به 
حجم زیادی از ماده خوراکی جهت آزمایش دارد. لذا امروزه استفاده از 
روش‌های برون‌تنی" توسط اکثر محققین ترجیح داده می‌شوند. اولین 
مطالعه در خصوص روش‌های برون تنی. مربوط به هضم دو مرحله 
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ای بود (۲۵). پس از آن. روش آزمون تولید گاز به منظور بررسی 
قابلیت هضم و محتوای انرژی مواد خوراکی مورد استفاده قرار گرفت 
(۱۴). در این روش مشاهده شد خوراک‌های مختلفی که در مایع 
شکمبه انکوبه شده بودنده مقدار گاز متفاوتی تولید کردند و بین حجم 
گاز تولید شده و ارزش انرژی‌زایی مواد خوراکی رابطه مستقیمی وجود 
داشت و این روش را آزمون تولید گاز" نام گذاری کردند (۲۳). روش 
آزمون تولید گاز اطلاعات با ارزشی درخصوص کینتیک هضم خوراک 
در شکمبه ارایه می‌دهد و با توجه به این که نیاز به تجهیزات ارزان 
قیمتی دارد لذا امروزه به‌عنوان یک تکنیک آزمایشگاهی, به‌طور 
وسیعی مورد پذیرش قرار گرفته است (۲۸). در روش آزمون تولید گازه 
حجم گاز تولید شده در زمان‌های متوالی مورد اندازه‌گیری قرار 
می‌گیرد و روند تولید گاز در آن به‌صورت یک منحنی غیر خطی ارایه 
می‌گردد. در یک منحنی تولید گاز ۳ بخش شامل فاز اولیه (فاز بدون 
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تولید گاز و یا تولید آهسته گاز)» فاز تولید گاز سریع و فاز پایانی (که در 
آن سرعت تولید گاز کاهش يافته و در نهایت به صفر می‌رسد) قابل 
تشخیص است. برای برآورد نتایج چنین کینتیکی نیاز به یک مدل 
ریاضی خاصی می‌باشد که بتواند روند تولید گاز را در طول زمان 
توصیف کند (۲). به‌عبارت دیگر توصیف ریاضی کینتیک تولید گاز 
(حجم گاز تولید شده در طول زمان) که از طریق برازش آن با تعدادی 
از مدل‌های ریاضی صورت می‌گیرد» اطلاعات با ارزشی در خصوص 
کینتیک تخمیر شکمبه ارایه می‌دهد و برای این منظور در طول 
دهه‌های گذشته تعداد زیادی از مدل‌های ریاضی توسط محققین 
مختلف پیشنهاد شده است (۷). مدل غیر خطی ریاضی که معمولا 
برای تعیین کینتیک تخمیر شکمبه استفاده می‌شود از نوع مدل 
نمایی" (۳26۳) می‌باشد که اولین بار در سال ۱۹۷۹ توسط اورسکوف 
رايج‌ترین مدل ها برای برازش داده‌های آزمون تولید گاز و پیش‌بینی 
کینتیک تخمیر شکمبه می‌باشد. اما این مدل ساختار ساده نیمه- 
هذلولی و غیر سیگوییدی دارد و چون فعالیت میکروبی در شکمبه و 
روند رشد آن‌ها به صورت سیگموییدی است. بنابراین استفاده از یک 
مدل غیر سیگموییدی مانند ۳26۳ برای تفسیر نتایج مربوطه. از دقت 
کافی برخوردار نیست و سبب ایجاد خطا در آنالیز داده‌ها می‌شود (۳). 
لذا» امروزه انواع مدل‌های با ساختار سیگموییدی و غیر سیگموییدی 
برای برآورد بهتر نتایج آزمون تولید گاز به کار برده شده‌اند و در این 
رابطه انواع معادلات مختلف ارایه شده است که هدف اصلی استفاده از 
این معادلات» پیش‌بینی دقیق‌تر کینتیک تخمیر شکمبه می‌باشد (۲۰). 
لازم به ذکر است که برخی از محققين هنوز هم ترجیح می‌دهند که از 
مدل ۲2۳ استفاده کنند. اما مدل‌هایی هم وجود دارند که برای این 
منظور ارایه شده‌اند اما کمتر مورد توجه قرار گرفته‌اند (۲۴). بر این 
اساس برخی از محققین استفاده از ساير مدل‌ها را ترجیح می‌دهند 
(۱۱). مدل لجستیک" (106) و مدل گومپرتز" (001) از پهترین 
برای توصیف کینتیک تولید گاز توسط محققین استفاده شدند (۲۲). 
مدل فیزو" (۳21) نیز به منظور بررسی کینتیک تخمیر شکمبه 
استفاده شد (۱۸). از آنجا که هر یک از مدل‌های پیشنهاد شده دارای 
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کارایی و قابلیت متفاوتی در مقایسه با یکدیگر می‌باشند. بنابراین آنالیز 
داده‌های به‌دست آمده از آزمایش آزمون تولید گاز به وسیله مدل‌های 
مختلف» ممکن است نتایج متفاوتی ارایه دهد. بر این اساس, تعیین 
تفاوت موجود در بین این مدل‌ها از نظر توانایی آن‌ها در پیش‌بینی 
بهتر نتایج» می‌تواند در انتخاب مدل ایده آل بسیار تأثیر گذار باشد (۱۶). 
برای مقایسه قابلیت مدل‌ها (نکویی‌برازش آن‌ها) در پیش‌بینی تتایج 
آزمون تولید گاز از ابزارهای مختلفی استفاده می‌شود. برخی از این 
ابزارهاء آماره‌هایی هستند که از تجزیه واریانس داده‌ها به‌دست 
می‌آیند. مانند میانگین مربعات خطاه (۷۲۳» ضریب تعیین * (2) 
انحراف مطلق میانگین باقیمانده" (0۷۸1) و میانگین درصد خطا؛ 
(01۳5). اما برخی دیگر از آماره‌ها از طریق آزمون‌هایی از قبیل آزمون 
دوربین -واتسون (650) طمداج۷۷-عذ0) و آزمون اجرا (50عا صنت) و 
نیز محاسبه رابطه رگرسیون خطی بین مقادیر مشاهده شده و با مقادیر 
پیش‌بینی شده محاسبه می‌شوند (۲۷). 

سیلاژ ذرت به‌عنوان یک خوراک علوفه‌ای رایج در تغذیه گاوهای 
شیری و پرواری می‌باشد. آگاهی از کینتیک تخمیر شکمبه‌ای این 
خوراک به‌خصوص در زمان‌های مختلف پس از سیلو کردن» سبب 
می‌شود که شناخت بهتری از این خوراک داشته باشیم. بنابراین هدف 
پژوهش حاضر ارزیابی نکویی برازش مدل‌های نمایی (۲۶۳) 
لجستیک (100) گومپرتز (00(4) و فیزو (211) به‌وسیله تکنیک 
آزمایشگاهی آزمون تولید گاز با سیلاژ ذرت به عنوان نمونه خوراک به 
کمک آماره‌های میانگین مربعات خطا ۷1818 ضریب تعیین (18 
انحراف مطلق میانگین باقیمانده (۷1۸) و ميانگین درصد خطا 
(۳۲) بود. 


مواد و روش‌ها 

آماده سازی نمونه‌های خوراک: در این آزمایش از سیلاژ 
ذرت به‌عنوان نمونه خوراک استفاده شد. این خوراک‌ها در زمان‌های 
صفر ۳۰ و ۶۰ روز پس از سیلو کردن تهیه شدند. ابتدا درصد ماده 
خشک نمونه‌ها با استفاده از آون در دمای ۶۰ درجه سلسیوس به‌مدت 
۸ ساعت تعیین شد. سپس نمونه‌های خشک شده با استفاده از 
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آسیاب رومیزی مجهز به الک ۱ میلی‌متری آسیاب شدند و از نمونه- 
های آسیاب شده برای تعیین ترکیب شیمیایی و انجام آزمون تولید 
گاز استفاده شد. ترکیب ۷ ۲ شیمیایی نمونه‌ها با استفاده از روش‌های 


استاندارد تعیین شد (۱ و ۲۶). نتایج مربوط به ترکیب شیمیایی سیلاژ 


ذرت در جدول ۱ ارایه شده است. 


جدول ۱- ترکیبات شیمیایی (برحسب درصد ماده خشک) سبلاژ ذرت در زمان‌های مختلف سیلو کردن 
ععصصتا عصتاتعصه اصمع1ل اه معملاه رمع ۶ (۵ 1۱۷۲ ۵ 860عه) ممت)لوممصورمی [همتصصمطزن -1 م12۳0 


لیاف تامتحلول ذر شوینده اسیدی. لیاف نامتحلول در شوینده ختفی 
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۲- آزمایش آزمون تولید گاز: آزمایش آزمون تولید گاز در ۲ 
اجرای (۲0) جداگانه و در طی ۲ دوره ۶ روزه انجام شد. برای اندازه- 
گیری تولید گاز از روش میزان جابجایی پیستون در داخل سرنگ 
شیشه ای مدرج مرتبط با ویال‌های ۵۰ میلی‌لیتری حاوی مایع شکمبه 
و نمونه خوراک. که معرف میزان تولید گاز بوده استفاده شد (۵). در 
هر اجرا» نمونه‌های خوراک در ۴ تکرار به همراه ۲ عدد ویال به‌عنوان 
بلانک در نظر گرفته شد. برای این منظور ابتدا مایم شکمبه از تعداد 
۳ رأس گوسفند نر مهربان مجهز به فیستولای شسکمبه‌ای قبل از 
خوراک‌دهی صبحگاهی گرفته شد. این گوسفندان روزانه در دو وعده 
غذایی صبح و عصر با استفاده از ۶۰ درصد علوفه (یونجه و سیلاژ 
ذرت) و ۳۰ درصد کنسانتره تغذیه می‌شدند (۱۵). مایع شکمبه به‌دست 
آمده از هر ۳ رأس گوسفند با هم مخلوط شدند و بلافاصله در شرایط 
بی‌هوازی و با استفاده از فلاسک در مجاورت آب گرم ۳۹ درجه 
سلسیوس به آزمایشگاه منتقل شد (۶). مایع شکمبه با استفاده از 
پارچه متقال چهار لایه صاف شد و به نسبت ۱ به ۲ با محلول بافر 
آماده شده در مجاورت گاز دیا کسیدکربن مخلوط شد و به‌عنوان مایع 
شکمبه بافری شده جهت استفاده بعدی در داخل بن‌ماری ۲٩‏ درجه 
سلسیوس نگهداری شد (۱۳). به‌داخل ویال‌های شیشه‌ای مربوط به 
نمونه‌های خوراک مقدار ۲۰۰ میلی‌گرم نمونه خوراک (خشک شده در 
آون ۶۰ درجه سلسیوس) به‌همراه ۳۰ میلی‌لیتر از مایع شکمیه یافری 
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کردن نمونه خوراک, فقط از مقدار ۳۰ میل‌لیتر مایع شکمبه بافری 
شده استفاده شد. پس از آماده کردن همه ویال‌های مربوط به هر اجرا؛ 
درب ویال‌ها با دقت با درپوش آلومینیومی پرس شد و بلافاصله به- 
داخل بن‌ماری ۲۹ درجه سلسیوس منتقل شد و حجم گاز حجم گاز 
تولیل فده زیر یی ی لیر در ستاغاف ‏ ان غا ۳ 
۴ ۷۲ ۰۹۶ ۱۲۰ و ۱۴۴ پس از انکوباسیون. با استفاده از یک 
سر سوزن متصل به یک سرنگ شیشه‌ای مدرج ثبت شد. 

۳- مدل‌های غیر خطی استفاده شده و نحوه مقایسه 
نکویی برازش آن‌ها: به‌منظور پیش بینی حجم گاز تولید شده در 
زمان‌های مختلف انکوباسیون از ۴ مدل غیر خطی استفاده شد. معادله 
مدل‌های مورد مطالعه در جدول ۲ ارایه شده است. برای این منظور» 
تعداد ۲۷ سری داده به‌دست آمده از آزمایش (۳ نوع نمونه خوراک» ۳ 
تکرار و در ۲ اجرای جداگانه) که مربوط به حجم گاز تولید شده در 
زمان‌های مختلف انکوباسیون بودند. استفاده شد. همه داده‌ها با 
استفاده از رویه 10وععتعع۲ عمعصنامها در نرم افزار 16.0 ٩۳98‏ 
5 107 به مدل‌های مورد نظر برازش شدند و نتایج به‌صورت 
میلی‌لیتر به‌ازای ۲۰۰ میلی‌گرم ماده خشک بیان شد. 

به‌منظور بررسی نکوبی‌برازش مدل‌ها از مقادیر میانگین مربعات 
خطا (۷15۳)», ضریب تعیین (157)» انحراف مطلق میانگین باقیمانده 
(۷۵۲) و میانگین درصد خطا (۷1۳) به‌دست آمده از هر مدل 
استفاده شد. 


۰ نشربه پژوهشهای علوم دامی ایران. جلد ٩۲‏ شماره۴ . زمستان ۱۳۹۹ 


جدول ۲- توصیف مدل های ریاضی غیر خطی استفاده شده در این مطالعه 
۲ ونطا صا موی کاعل مجط 1هعن)معمعط/هدط مممصتلطمط ماه ت4صهه ۵۶ ممتاوتن وع1 -2 ماماه 1 


دامنه پارامتر ساختاری معادله مدل‌ها 
صتقصو(1[ ۲عاعصصه‌تهم ممرف٩‏ ۲ ممتامتاوظ۴ 015( 
0 ِ (*۳ ۰ -۸۰)1 < ۲ مدل نمایی 
(۸۳ظ) اهتادوممومد۴ 

۱20 9 رو ۵,1 -< ۲ مدل فیزو 
(۳2۲۷) طفیاطماز۲ 

20 0 ۶8 اي .یی ع و مدل گومپرتز 
(00(۳۷) 2اوحجط۵ت) 

120 9 دی + 1) .۵ < ۲ مدل لجستیک 
(1۷00) عتاعنعم] 


۱ ۷ : حجم گاز تولید شده در زمان ]» ۸ : پتانسیل تولید گاز 6 : سرعت تولید گازه 0: پارامتر ساختاری» : زمان انکوباسیون و 6 : عدد نیر (.۲/۷۱۸۲۱۸۲۸۸۴) 


احماعجمی و توملا :ع 4مه عصصتا حمتاعطبمصا ۰ رتماعمتهم ممقطه :0 رتمتمصمتمم ملق نع رمصامه قمع متاماموطرروه ۸ با متا کم قمع ۵۴ میامن 1 


مقادیر ۷5۳ و 2 در هر مدل از رابطه‌های ۰۱ ۲ و ۲ محاسبه 


شدند (۵): 
رابطه ۱ *( - 2 ع ۲25 
ت 
- 158 
زانط ۷ :۲ 
۱ 
--1-۳- ۲3 
رابطه ۳ 5 


در رابطه‌های ۱ ۲ و ۰۳ ۷ مقدار مشاهده شده» # مقدار پیش- 
بینی شده 165 محموع مربعات باقیمانده ظ تعداد داده‌ها (تعداد نقاط 
زمانی اندازه‌گیری شده)» 0 تعداد پارامتر موجود در هر مدل و 517 به- 
عنوان واریانس کل در مقادیر مشاهده شده می‌باشد. 

مقادیر ۷۸ و ۷۳۳۲ نیز بر اساس رابطه‌های ۴ و ۵ به‌صورت 


ار تزا 


۳ 
رابطه ۴ 8 وسزط < ۲۳۸۲ 
من تاج 2 مور 
رابطه ۵ [ ۲ 


بررسی قابلیت مدل‌ها در برازش داده‌ها با استفاده از آزمون 
دوربین -واتسون و آزمون اجرا (1651 عد) استفاده شد (۴ و ۲۷). رابطه 
رگرسیون خطی بین مقادیر مشاهده شده (محور ) و پیش‌بینی شده 
(محور ۷) در هر یک از زمان‌های انکوباسیون برای همه سری داده‌ها 
و در هر مدل, بر اساس رویه 1۳0 با استفاده از برنامه ٩۸٩‏ (۲۱) 
محاسبه شد و مقادیر عرض از مبداء (101676608) و ضریب رگرسیون 


(51026) برای هر مدل به‌دست آمد. معنی‌دار شدن پارامترهای 


)2.718218284...( 


رگرسیونی نیز بر اساس فرضیه 10100601-0 و ٩10061‏ در سطح 
خطای ۵ درصد آنالیز آماری شدند (۲۷). 

احتمال معنی‌دار بودن تفاوت مقدار 101670601 به دست آمده از 
مدل‌ها با عدد صفر از طریق رویه ۳0 و با استفاده از برنامه 5۸6 
(1999) انجام شد. احتمال معنی دار بودن تفاوت مقدار 51006 با عدد 
۱ نیز از طریق آزمون ) مطابق رابطه ۶ محاسبه شد و مقدار ] به‌دست 
آمده از هر مدل» با ] جدول و در سطح خطای ۵ درصد مقایسه آماری 
شد : 


1 


25 5 ت 


رابطه ۶ نود 


در رابطه ۶ «ا مقدار ) محاسبه شده برای هر مدل, :0 مقدار 
ضریب رگرسیون محاسبه شده در مدل ء 1/1515 میانگین مربعات خطا 
و 526 مجموع مربعات مقادیر مشاهده شده می‌باشد. 

- آنالیز آماری داده ها: آنالیز آماری مقادیر برآورد شده از 
طریق هر کدام از مدل‌های مورد مطالعه (شامل پارامترهای ۸۵ و ع) و 
پارامترهای مربوط به نکویی‌برازش (۷151, 67 14۸ و ۳۳1) 
بر اساس رویه 01/۷ با استفاده از برنامه ٩۸٩‏ (۲۱) انجام شد. مدل 
آماری استفاده شده به‌صورت زه +13 1۷:۲ جر ع ۷ بود که در آن 
۷ متغییر وابسته. ‏ میانگین» :۷1 اثر مدل. 16 اثر ران و زه خطای 
باقیمانده بود. مقایسه میانگین‌ها نیز با استفاده از آزسون توکی و در 
سطح خطای ۵ درصد انجام شد. 


زابلی و همکاران,بررسی دقت برخی از مدل‌ها در برآورد ضرایب آزمون تولید گاز در سیلاژ ذرت ۳۷۱ 


جدول ۳- مقایسه فراسنجه‌های کینتیک تخمیر شکمبه‌ای برآورد شده بوسیله مدل‌های مورد مطالعه 
کام0مجط 0160و عط بوها 4متمصتاوی کم‌امصصه‌تومر متامصت! جمت)‌هاهم‌صررع؟ [مصرتصصیار ما ۵۶ طموزتهی‌مجمن -3 6اه 1 


مقدار ‏ مدل لجستیک مدل گومپرتز مدل فیزو 
7 ما۷۵ مناعلع 1[ ۹ طوباطماز۲ 
234«( 910 9837 914 
001> *05. 0 02053 064 09" 
۰ 1048 1465 2.1 
2.519«( 9127 9۹28 913 
0.001> 098۹ 0 32۰«( 304«( 
۰ 1164 1368 4 
906« 76.09 1123 165 
04 ۳ 0.093 0062 0.05 
. 1402 1512 2.12 
۹313 0 96952 1196 34 
4012 0 0۹5 0 59« 0.011 
. 1205 1615 ۳/۳ 


مدل نمایی پارامتر زمان بعد از سیلو کردن 
امتمعجويدظ .۰( ممعاعهمنوط عصنااممی ماه وعصز] 
9616 ۸ روز صفر 
۰057 0 ۰ 12۷0 
3 9 
119 ۸ روز ۳۰ 
رد0۳ ۰ 0 127 
9 9 
77.18 ۸ روز ۶۰ 
00115 6۵ 0 127 
۳ 9 
95994 ۸ میانگین 
0.017 ۰ ۷2( 
۴ 9 


" ۸: پتانسیل تولید گاز (میلی‌لیتر بر ۲۰۰ میلی‌گرم ماده خشک) «ا : پارامتر ساختاری» 06 : سرعت تولید گاز (میلی‌لیتر بر ساعت) 
۲ میانگین‌های دارای حروف متفاوت در هر ردیف از جدول دارای تفاوت معنی‌دار (0>۰/۰۵) بر اساس آزمون توکی هستند. 


6 متمصمتمم ملق جع رتعاعصمهم ممقطه :0 رر(ظ مه200/ل) عنام قمع متامام‌رروه خه 1 
۵۰ ۶ اوما عطا بو (۳>۵۰۵05) عصمتماگنل راامم‌تاوناهاه مه واونموتمونه کصعتعلکنک توا فععملام؟ مهوت منطانه مفعا۱ 2 


نتایج مربوط به پارامترهای پیش‌بینی شده توسط مدل‌های مورد 
مطالعه شامل پتانسیل تولید گاز (۸ پارمتر ساختاری" منحنی (0) و 
سرعت تولید گاز (0) در جدول ۳ ارایه شده است. مدل‌های مورد 
مطالعه از نظر پیش‌بینی مقدار پارامتر ۸ تفاوت معنی‌داری با هم نشان 
ندادند و مقدار پارامتر ۸ در بین مدل‌های مورد مطالعه در دامنه 
۷ (در مدل 003۷ در روز صفر) تا ۷۶/۰۹ (در مدل 10 در 
روز ۶۰) میلی‌لیتر به ازای ۲۰۰ میلی گرم ماده خشک قرار داشت. لازم 
به ذکر است که پژوهش مشابهی در خصوص استفاده از مدل‌های 
مورد بررسی در این مطالعه بر روی سیلاژ ذرت» در دسترس ما قرار 
نگرفت. اما در یک مطالعه که بر روی ۲۲ نوع خوراک علوفه‌ای 
مختلف انجام شده بود. مقدار پارامتر ۸ پیش‌بینی شده توسط 
مدل‌های ۰۳26۳ 00۷ و 100 در طول ۴۸ ساعت انکوباسیون» 
تفاوت معنی‌داری با هم نشان دادند (۲۷). احتمالاً ترکیب شیمیانی 
خوراک‌های مورد آزمایش شرایط انجام آزمایش, جمعیت میکروبی 
مایع شکمبه و زمان انکوباسیون از مهم‌ترین عوامل تأثیر گذار در اين 
خصوص می‌باشند. 

در تعیین ارزش غذایی مواد خوراکی در آزمون تولید گاز » حجم 
گاز تولید شده پایانی و پتانسیل تولید کاز بسیار مهم است. پتانسیل 
تولید گاز در مواد خوراکی به فاکتورهای متعددی از قبیل مقدار نمونه. 


اهتدم م2 - [ 


اندازه ذرات» دوره انکوباسیون, نوع جمعیت میکروبی مایع شکمبه و 
نیز شرایط آماده سازی مایع شکمبه دارد. لذا انتظار تولید مقادیر مشابه 
گاز تولیدی در زمان‌های مختلف انکوباسیون و در آزمایش‌های 
مختلف امکان پذیر نیست (۲۳). 

سرعت تولید گاز (0) در بین مدل‌های مورد بررسی فقط در 
زمان‌های صفر و ۲۰ سیلو کردن, دارای تفاوت معنی‌دار بود و مقدار 
آن در مدل 1006 بیشتر از سایر مدل‌ها بود (۳>۰/۰۵). اما مقدار 
پارامتر » در میانگین کل دوره سیلو کردن در بین مدل‌ها تفاوت 
معنی‌داری نشان نداد. در مطالعه‌ای که بر روی سیلاژ ذرت انجام 
گرفته بود مقدار عددی ‏ برآورد شده از طریق مدل ۳2۳ مقدار 
۵ میلی‌لیتر بر ساعت گزارش شد (۸). همچنین در مطالعه دیگره 
مقدار آن ۰/۰۴۳ میلی‌لیتر بر ساعت گزارش شد .)۱٩(‏ 

با توجه به اينکه پارامتر ۵ یک پارامتر ساختاری است که برای 
تصحیح منحنی به مدل اضافه می‌شود و از آنجا که مقدار آن با توجه 
به ساختار معادله در مدل‌های مختلف تغییر می‌کنده لذا امکان مقایسه 
کردن پارامتر ‏ در بین مدل‌های مختلف وجود ندارد (۲۹). در بین 
مدل‌های مورد مطالعه» مدل 126۳ یک مدل ساده نمایی است که 
فقط از دو پارامتر ۸ و » تشکیل شده است و فاقد پارامتر ط می‌باشد. 
اما سایر مدل‌ها دارای پارامتر 0 می‌باشند. وجود پارامتر 9 در سایر 
مدل‌ها در مقایسه با مدل 726۳ احتمالاً سیب افزایش انعطاف پذیری 


آن‌ها در برازش داده‌ها و تصحیح روند تولید گاز شده و در نهایت 


۲۳ نشربه پژوهشهای علوم دامی ایران. جلد ٩۲‏ شماره۴ ۰ زمستان ۱۳۹۹ 


۲ در مدل ۳6۳ و کمترین آن در مدل 001۷ مشاهده شد. 
مقدار 122 نیز در مدل‌های 122۳ 001 و 106 در آن مطالعه به- 
ترتیب ۰۹۶۹ ۰/۹۹۲ و ۰/۹۸۵ به‌دست آمد (۲۷). در مطالعه دیگری 


باعث افزايش دقت آن مدل‌ها گردیده است (۲۷). 


نتایج مربوط به پارامترهای نکویی برازش در بین مدل‌های مورد 
مطالعه در جدول ۴ ارایه شده است. مطابق جدول فوق. مدل ۳۲۳ 


دارای بیشترین مقدار ۷۹۳ و کمترین مقدار 27 در بین مدل‌ها بود. که بر روی ۱۵ نوع سیلاژ خانواده گندمیان به روش آزمون تولید گاز 


اما بر خلاف آن» مدل 003۷ و ۳211 دارای کمترین مقدار ۷۹۳ و 


در مدت زمان ۷۲ ساعت انکوباسیون صورت گرفت. مقدار ۷18۳۲ در 
بیشترین مقدار 1*2 بودند. در یک آزمایش مشابه که بر روی مواد مدل‌های 1226۳ 003۷ و 100 به ترتیب ۰۲۳/۸ ۱۰/۲ و ۲۵/۸ و 
خوراکی مختلف به‌روش آزمون تولید گاز انجام شد بالاترین مقدار مقدار 2 نیز به ترتیب ۰۰/۹۹۸ ۰/۹۹۹ و ۰/۹۹۷ گزارش شد (۷). 


جدول 6- مقایسه مدل‌ها بر اساس پرامترهای نکوبی برازش آنها 
منهج 1 01 55عظ000م تتعط) چرم 00وعه ورم0ممطر مط 0۶ طموزتجمرمومن -4 ماماه 1 


مقدار و مدل لجستیک مدل گومپرتز مدل فیزو مدل نمایی پارامترهای نکوبی برازش 
مبااه و 100 اتازفرگ ۳ 5( مهن 1 اه موعملوون0 
01> 25.4 12 2-06۳ 292-30 میانگین مربعات خطا 
2 9۳( 

01> 90 0۳ 990 0 ۳ 0975 00961 ضریب تعیین 
# 2 

001),> 2-2-153۰ ۴ 2-91 29۳ 22-2073 انحراف مطلق میانگین باقیمانده 
۳/۸۵ 

001),> ۲ ۳ 5 ۳ میانگین درصد خطا 
۷۳۴۶ 


" میانگین‌های دارای حروف متفاوت در هر ردیف دارای تفاوت معنی‌دار (0>۰/۰۵) بر اساس آزمون توکی هستند. 


۵۸ 0 )وه مطا بط (۳>۵۰۵05) اصمتمزگنل برالمع‌تامتاهاه میج واوتموتمم‌نه کصمتقتگنل رها 0ععملام؟ مسمد صنطانه ممع۱ ۱ 


منحنی تولید گاز یک شکل لگاربتمی دارد و به همین دلیل مدل- 
هایی که برای توصیف آن به کار می رونده اکثرً فرم لگاریتمی و غیر 
خطی دارند. مقایسه مدل‌های غیر خطی از طریق تکنیک نکویی 
برازش و آنالیز باقیمانده (مقدار مشاهده شده منهای مقدار پیش بینی 
شده) يا همان مقدار خطا صورت می‌گیرد (۱۲). بر اساس نظر 
محققین, مقدار 12 در مدل‌های توصیف کننده آزمون تولید گاز معیار 
مناسبی برای قضاوت و انتخاب بین مدل‌هاست (۲۹). 

گزارش گردیده است که مدل‌های با تعداد پارامتر بیشتر معمولا 
مقدار 1۷1512 کوچکتری دارند که در نتایج ما نیز چنین روندی مشاهده 
شد (۲۰). گزارش شده است که اگر مقدار 1۷۲۳ و 7 در یک مدل 
به‌ترتیب زیاد و کم باشد. نشان دهنده این است که مدل دارای 
نکویی‌برازش خوبی نیست (۱۰ و ۱۲). بنابراین با توجه به 
مقادیر ۹۳ و 32 به دست آمده» مشخص شد که مدل ۳۳ 
به‌خاطر بیشسترین مقدار 151" و کمترین مقدار 18 ضعیف‌ترین 
عملکرد را در بین همه مدل‌ها داشت. 

مقدار 1۸ نیز به‌عنوان یک شاخص برای بررسی نکویی- 


کتمته ومتمناوو صقمصر :۱6۴ 2 
ی 
صمناهع عانامعطاه صفعجه آمبل‌تفور :۸۲( 3 
عم مهمامه‌هم صعمصط :۱۳۴۲ 4 


برازش مدل‌ها استفاده می‌شود. هر چقدر مقدار 5۷۲۸۲ کوچکتر 
باشد (به صفر نزدیک‌تر باشد) نشان‌دهنده قابلیت بهتر مدل در برازش 
داده‌هاست (۱۶). مطابق جدول ۴ مقدار 1/۸ در مدل‌های 
و ۳2۳۲ کمترین مقدار (به‌ترتیب ۲/۵۹۱ و ۲/۸۷۹) و در مدل 
۳ بیشترین مقدار (۲/۸۰۷) بود (0>۰/۰۵). بر این اساس» 
مدل‌های 00۷۲ و ۳2۲ در مقایسه با مدل ۳26۳ قابلیت بهتری در 
برازش داده‌ها داشتند. 

مقدار ۷1۳5 نیز نشان‌دهنده این است که آیا میانگین مقادیر 
پیش‌بینی شده نسبت به مقادیر مشاهده شده. بیش برآورد و یا کم 
برآورد شده است یا خیر. بر اساس پارامتر ۷1۳12 اگر مقادیر پیش‌بینی 
شده توسط مدل نسبت به مقادیر مشاهده شده بیشتر باشد (بیش 
برآورد باشد) مقدار ۷1۳[ منفی و در غیر این صورت مثبت است. 
همچنین. اگر مقدار این عدد به صفر نزدیک‌تر باشد. نشان‌دهنده بهتر 
بودن نکویی‌برازش مدل می‌باشد (۱۶). در جدول ۴ مقدار ۱۷۳1 در 
مدل‌های 126۳ ۲211 00 و 00 به‌ترتیب ۰۵/۵۲۷ ۲/۰۷۵ 
۸- و ۳/۸۲۰- بود و مقدار آن در مدل ۳226۳ تفاوت معنی‌داری 


زابلی و همکاران,بررسی دقت برخی از مدل‌ها در برآورد ضرایب آزمون تولید گاز در سیلاژ ذرت ۳۷۳ 


با سایر مدل‌ها نشان داد (۳>۰/۰۵). بر این اساسء مقادیر پیش بینی 
شده توسط مدل ۲2۲۳ کمتر از مقادیر مشاهده شده بود (کم براورد 
شده است). اما سایر مدل‌ها هر چند که مقدار اندکی بیش‌برآورد 
داشتند» اما وضعیت آنها در این رابطه بهتر از مدل ۳726۳ بود. 

در رابطه با مقادیر 31۸ و ۷۳ به‌دست آمده از مدل‌های 
مورد مطالعه در سیلاژ ذرت» گزارشی در دسترس ما قرار نگرفت. اما 
مشابه نتایج ماء در مطالعه‌ای که بر روی علوفه پونجه صورت گرفت؛ 
مقدار 1371۸ و ۷۳۲ در مدل ۳26۳ به‌ترتیب ۲/۲۶۲ و ۲/۷۸۶ به- 
دست آمد که تفاوت معنی‌داری با سایر مدل‌ها داشت. این وضعیت 
نشان‌دهنده قابلیت ضعیف‌تر مدل ۳26۳ در برازش داده‌های گاز تست 
برای علوفه یونجه بود (۲۹). همچنین» در مطالعه‌ای که بر روی یک 
جیره حاوی ۶۰ درصد یونجه و با مدل‌های مختلف انجام شد مقدار 
عددی ۸۲ و ۷۳ در مدل 1726۳ بیشتر از مدل‌های دیگر بود 
که نشان‌دهنده نکویی‌برازش ضعیف‌تر مدل ۳ در مقایسه با 
مدل‌های دیگر بود (۱۶). 

نتایج مربوط به آزمون دوربین-واتسون (99) در جدول ۵ ارایه 
شده است. دامنه اماره 1۷۷ در محدوده عددی ۰ تا ۴ می‌باشد و هر 
چقدر مقدار 1۷۷ به‌دست آمده از یک مدل به مرکز این محدوده 
(عدد ۲) نزدیک‌تر باشد. نشان‌دهنده پراکندگی مستقل خطاها (مقادیر 
باقیمانده) می‌باشد و در غیر این صورت و در اعداد نزدیک به مرز این 
محدوده ( اعداد صفر و ۴) خطاهای به دست امده از مدل دارای خود 
همبستگی هستند (۲۹). مطابق جدول ۵ در تمامی مدل های مورد 


مطالعه» مقدار آماره 12۷7 به عدد صفر نزدیک بود. این وضعیت در 
مدل ۳26۳ ملموس تر بود که نشان داد ضعف مدل 726۳ در برازش 
داده‌ها بیشتر از سایر مدل‌ها است. همچنین در تمامی مدل‌های مورد 
بررسی» فرض 110 رد شد و با توجه به اینکه رد شدن فرض ۳ 
نشان‌دهنده وجود خود همبسنگی بین خطاها (مقادیر باقیمانده) می- 
باشد. لذا وجود خود همبستگی در بین خطاها از این طریق نیر تأیید 
شد و نشان داد که همه مدل‌ها در برازش داده‌ها دارای ضعف بودند. 
اما ضعف مدل 720۳ در اين رابطه بیشتر بود. بر اساس نظر محققین» 
معنی‌دار شدن آماره 1۷۷ نیز نشان دهنده وجود خود همبستگی در 
بین خطاها و ضعف مدل در برازش داده‌ها می‌باشد (۲۹). همانطور که 
در جدول ۵ مشاهده می‌شود. از بین تعداد ۲۷ منحنی بررسی شده 
برای هر مدل, بیشترین تعداد معنی‌دار شدن آماره ۷۷ در مدل 
6۳ (در ۲۴ منحنی) و کمترین آن در مدل 00۷ (در ۱۸ منحنی) 
مشاهده شد که تأیید کننده نتایج فوق می‌باشد. در مطالعه‌ای که بر 
روی ۲۳ نوع اقلام خوراکی مختلف در قالب ۲۳ منحنی تولید گاز و با 
استفاده از مدل‌های ۳26۳ 00 و 1,0 صورت گرفت» مشخص 
شد که در مدل 326۳ در تعداد بیشتری از منحنی‌ها آماره 1۷۷ 
معنی‌دار شد (تعداد ۲۲ منحنی). اما مدل 00۷1 کمترین مقدار (تعداد 
۶ منحنی) را در اين رابطه داشت (۲۷). لازم به ذکر است که وجود 
خود همبستگی بین مقادیر خطاها در یک مدل به تنهایی نمی‌تواند 
ملاک ارزشیابی یک مدل از نظر نکویی‌برازش آن باشد و نیاز است تا 
سایر شاخص‌ها نیز بررسی گردد (۱۷ و ۲۷). 


جدول ۵- ننایج آزمون دوربین-واتسون در مدل‌های مورد مطالعه 
اواعمع 0عن0ج0ه عط) صا واانافع؟ وع] (۲۱۱۷) جوفنه ۷۷ صتطظ -5 عاطو] 


مدل لجستیک مدل گومپرتز مدل فیزو 
۱996 اتازفزگی ۳71 
0۳615 0۳05 0691 
۹۵0 ۹۵0 ۹0 


100 


مدل نمایی 
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0392 آماره و رش وا تون 


۱7 
نتیجه فرض 0 
" عمط ما[ 


توزیع منحنی ها بر اساس معنی‌دار بودن و يا معنی‌دار نبودن (0<۰/۰۵) مقادیر آماره دوربین-واتسون (در همه ۲۷ منحنی) 
هام ومبطاه۷ ۲۱۷۷ (0.05<م) صهم)نصع1و-دمه 4صح (0.05>م) اصه‌تگتصهنه مطا مت مصت 26601 وعسن مط ۵۶ ممتاتاماتتاوزر[ 


(27عو07ه 
23 18 21 24 معنی دار شده 
اصهم1 نم وزو 
0 0 0 0 معنی دار نشده 
۱۹ 
4 9 6 3 بی معنی 


۲ فرض صفر: هیچ خود همبستگی معنی‌داری بین خطاها (مقادیر باقیمانده) وجود ندارد. 


معط[ 


۰(ومیاه۷ آمیل‌توع) ومع معط عصمصصه صمتاهاهتمممفنيه اصمعز/تصع زو ه امط و 16۲6 :و 1 


۴ نشریه پژوهشهای علوم دامی ایران. جلد ۲ شماره۴ ۰ زمستان ۱۳۹۹ 


در تحقیق حاضر نیز هر چند در مدل‌های 001۷ و 100 در 
بین خطاها خود همبستگی مشاهده گردید اما مطابق نتایج مربوط به 
نکویی‌برازش, این دو مدل قابلیت خوبی در برازش داده‌ها در مقایسه 
با مدل‌های 1226۳ و ۳2۳۷ داشتند. 

نتایج مربوط به آزمون اجرا (1691 عات) در جدول ۶ ارایه شده 
است. این آزمون به‌منظور بررسی روند تغییر داده‌ها (مقادیر باقیمانده) 
در هر مدل استفاده می‌شود و تغبیرات سیستمیک مقادیر فوق به‌طور 
واضح در نتایج آزمون اجرا مشاهده می‌شود (۳). یک اجرا (صنات)» تعداد 
اعداد پشت سر هم با علامت مشابه (مثلا + یا -) در یک سری اعداد 
متوالی می‌باشد. بر این اساس, اگر تعداد اجراهای موجود در مقادیر 
باقیمانده در یک مدل بیشتر باشذء احتمالا تکویی‌برازش آن مدل 
بیشتر است (۲۷). در مطالعه ماء مطابق جدول ۶ تعداد اجراهای به- 
دست آمده از ۲۷ منحنی که برای برازش هر مدل استفاده شده بودند. 
در ۴ دسته تقسیم شدند. این دسته‌ها شامل تعداد اجراهای با تعداد ۳ 
> ۰۴-۵ ۶-۷ و ۸< اجرا بودند. وقتی که مقادیر باقیمانده از هر 
مدل به‌طور تصادفی پراکنده باشند» تعداد اجرای مشاهده شده بیشتر 
خواهد بود و مدل نکویی برازش قوی‌تر و خطای سیستمیک کمتری 


خواهد داشت .)٩(‏ مطابق جدول فوق. در مدل 1726۳ همه ۲۷ سری 
داده‌ها (منحنی‌ها) دارای کمترین اجرا (۲ ک) بودند که نشان‌دهنده 
عملکرد ضعیف این مدل بود. سایر مدل‌ها به مراتب عملکرد بهتری 
نسبت به مدل 1726۳ داشتند و بیشتر منحنی‌های آنها دارای اجرای 
۴-۵ بودند. مشابه نتایج ماء در مطالعه‌ای که بر روی اقلام مختلف 
مواد خوراکی صورت گرفت» مدل ۳2۳ در مقایسه با مدل‌های 
7 و 0 از نظر آزمون اجرا عملکرد ضعیف‌تری داشت و مدل 
7 دارای بهترین عملکرد بود (۲۷). همچنین در مطالعه دیگری 
که بر روی اقلام خوراکی مختلف از جمله سیلاژ ذرت انجام شد از 
مجموع ۲۱۶ منحنی مطالعه شده. تعداد ٩۰‏ درصد و يا بیشتر منحنی - 
ها در مدل‌های ۰۳26۳ 003 و 100 دارای تعداد اجرای ۵ > بودند 
(۲). گزارش گردیده است که تعداد اجرای کمتر در یک مدل نشان- 
دهنده آن است که در مدل بیش‌برآورد و یا کم‌برآورد اتفاق افتاده 
است. در حالی که تعداد اجرای بیشتر نشان دهنده پراکنش تصادفی 
مقادیر باقیمانده (خطاها) بوده و بنابراین توانایی بهتر مدل ر توصیف 
می‌کند (۲). 


جدول - ننایج تست اجرا پس از برازش مدل‌ها 
داعم عط) عصتاز؟ ماه فالناعع اوعا سب -6 12016 


مدل نمایی 
۱۹9 


توزیع منحنی‌ها بر اساس تعداد اجراها (در همه ۲۷ منحنی) 
(۳7۵-27 آهاما) عصتا ۵۶ تعحاحصاط مطا مه عصتتمعمه ووسیه ۵۶ طمزاناماتتاوزنز 


مدل لحستیک مدل گومپرتز مدل فیزو 
100 /00 ۳11 
0 0 
26 25 26 
1 2 
0 0 


27 تعداد اجراهای کمتر از ۳ 
3> 
0 تعداد اجراهای بین ۴ تا ۵ 
4-5 
0 تعداد اجراهای بین ۶ تا ۷ 
67 
0 تعداد اجراهای بیشتر از ۸ 
8< 


توزیع منحنی‌ها بر اساس معنی‌دار بودن و معنی‌دار نبودن (0>۰/۰۵) نتایج تست اجرا (در همه ۲۷ منحنی) 
1) ات۲۵5۱ ]۵5] با ۵1 (0.05<) کصهعزگنطع1وحطمه 4صح (0.05>م) عصهم‌تگتصهنه مط ما مصتل01عع2 وعبسه مطا ۵۶ طمتاناماتتوزر[ 


راگیس هه هت ی یی قاس ماع خسف که سا یک امد تم شگی سسکا 
3 16 13 24 تفه 

" اصمهمز/نمعز۹ 

18 11 14 3 معتی دار تفنده 


۱ 
" فرض صفر: توزیع داده‌ها به صورت تصادفی است. 
اگر 2۱-۰2 < |2] باشه فرض صفر پذیرفته می‌شود و در غیر این صورت. فرض صفر رد می‌شود (227/۶۹۳۸ و ۰/۰۵ 0) 
حد بحرانی ,2 : ۱/۹۶ < 72۱-2 ناحیه بحرانی :اگر 1.96 < |2/ باشد» فرض صفر رد می‌شود. 
از آنجاییکه آماره تست (2) بیشتر از مقدار بحرانی است, نتیجه گیری می‌کنیم که داد ها در سطح خطای ۰/۰۵ به صورت تصادفی توزیع نشده اند. 

,0ص میج منم ۵۶ صمتنباطاتطون 16 : ول ! 
8 << 72 تادناهاه اوع] م ۲6[666 ,50زتعطاه مت 200001 روم2۱ < 2۱ 1۶ 

6 < 2۸۱ ]1 ۲0 ۲6661 :صمنعع آمعتاتن 6 ع 2۱ :(انها عوونا) فتاه آهعنانن 
۱6266 تمصع نی 0.05 عمط اج حمی‌هه۲ )مظ مج هاع مط) اقطه عمجم 9۵ رمباله۷ لمعناتت مطا صقطا تماهمته وا (2) متاوتاهاه اوعا عطا مممزگ 


512011168006 16761: 0 2 5 


زابلی و همکاران,بررسی دقت برخی از مدل‌ها در برآورد ضرایب آزمون تولید گاز در سیلاژ ذرت ۳۷۵ 


معنی‌دار شدن آزمون اجرا در یک مدل نشان‌دهنده آن است که 
ترتیب داده‌ها (مقادیر باقیمانده) دارای توزیع تصادفی نیستند و لذا 
نوعی خطای سیستمیک در مدل وجود دارد (۲۷). مطابق جدول فوق» 
از تعداد ۲۷ منحنی بررسی شده در مطالعه حاضر در مدل ۲2۳ در 
مقایسه با سایر مدل‌ها در تعداد بیشتری از منحنی‌ها (۲۴ منحنی) 
نتیجه آزمون اجرا معنی‌دار شد و مدل‌های ۳211 /00(۷ و 106 
وضعیت بهتری در مقایسه با مدل ۳26۳ داشتند که مطابق با نتایج 
سایر محققین می‌باشد (۲۷). 

نتایج مربوط به رابطه رگرسیونی بین مقدار حجم گاز مشاهده 
شده و مقدار پیش‌بینی شده آن در طول زمان انکوباسیون در جدول ۷ 
ارایه شده است. معنی‌دار شدن پارامترهای رگرسیونی بر اساس فرضیه 
۰< 101070601 (عرض از مبداء) و ۱< 91006 (ضریب رگرسیون) آنالیز 
آماری شدند. در چنین شرایطی بیشترین همبستگی بین مقدار 
مشاهده شده با مقدار پیش‌بینی شده در حالتی به‌دست می‌آید که 
مقدار عرض از مبداء و ضریب رگرسیون به ترتیب تفاوت آماری 
معنی‌داری با اعداد صفر و ۱ نداشته باشند (۱۷ و ۲۷). همان‌طور که 
در جدول ۷ مشاهده می‌شود. در همه مدل‌ها مقادیر عرض از مبداء و 
ضریب رگرسیونی به‌ترتیب تفاوت معنی‌داری با اعداد صفر و ۱ داشتند. 
مقدار عرض از مبداء در مدل ۳26۳ منفی و در سایر مدل‌ها مثبت بود. 


بنابراین مطابق شکل ۱ این وضعیت سبب شد تا خط رگرسیونی 
برآورد شده در مدل 1726۳ به خصوص در اوایل دوره انکوباسیون در 
زیر خط ۱:۱ قرار بگیرد. در چنین وضعیتی مقدار گاز پیش‌بینی شده 
(به خصوص در ساعات اولیه انکوباسیون) کمتر از مقدار مشاهده شده 
می‌باشد. این وضعیت از نظر بیولوژیکی امکان‌پذیر نیست. زیرا حجم 
کار یفنم در ساعت عفر تکریایون هقی رم شوه 32۳0 
سایر مدل‌ها (مدل‌های ۳217 001 و ۲00) عکس این حالت 
مشاهده شد و خط رگرسیونی برآورد شده در ساعات اولیه انکوباسیون 
در بالای خط ۱:۱ قرار گرفتند و سبب شد که حجم گاز پیش‌بینی 
شده در ساعت صفر انکوباسیون مثبت شود. از آنجا که در حالت 
ظبیعی بحتم کار تلی ده در متفر انکوباسیوی صقر من 
باشد. بنابراین بروز چنین وضعیتی نیز از نظر بیولوژیکی امک‌ان‌پذیر 
نیست و همه این‌ها دلالت بر عدم دقت زیاد (۱۰۰ درصد) در مدل- 
های مربوطه دارد. با این حال در بین مدل‌های مطالعه شده خط 
رگرسیونی برآورد شده توسط مدل‌های ۳211 و 003 به خط ۱:۱ 
نزدیک‌تر بود و مطابق جدول ۶ مقادیر :۷15۲ و *18 نیز در آن‌ها به- 
ترتیب کوچک‌تر و بزرگ‌تر از سایر مدل‌ها بود و بنابراین مدل‌های 
71 و 00 در مقایسه با سایر مدل‌های مورد مطالعه پیش‌بینی 


بهتری از پروفیل تولید گاز داشتند. 


جدول ۷- پارامترهای رگرسیونی بین حجم گاز مشاهده شده (در محور 26) و پیش بینی شده (در محور ۷) بعد از برازش مدل‌ها 
0۰ 6ظ) عصناز1 تماله ومصصتاو۷ عقع رفن2 ۷ 0 0ممن0عنرور فنایزع۷ (26-229 6 ۷۵0 عوها۵ جمعهاع وتعتعصصوتهن مماووع/ع1۵ -7 12016 


ضریب تعیین میانگین مربعات خطا ضریب رگرسیون 
٩06۲ (۲ ۳۳‏ 
0962 32006 1046 
00917 15-577 0969 
09۹1 120468 *093 
0912 1-9-۹414 0959 
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نتیجه گیری کلی 


تایج این آزمایش نشان داد مقدار پتانسیل تولید کاز (۸) پیش 
بینی شده توسط مدل‌های مورد مطالعه. تفاوت معنی‌داری با هم 


نداشتند. اما بررسی پارامترهای نکویی برازش نشان داد که در بین 
مدل‌های بررسی شده مدل 00۷ و ۲211 در مقایسه با مدل ۲2۳ 


کینتیک تخمیر شکمبه‌ای سیلاژ ذرت را با دقت بیشتری پیش‌بینی 
کردند. 


۴۷۶ نشریه پژوهشهای علوم دامی ایران. جلد ۱۲ شماره۴ ۰ زمستان ۱۳۹۹ 
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۰ رطهاعصتاعظ روافتصعطی آهه‌تاراممد لههه‌تلان ۵۶ ممتاهاه0وفظ .60 طا16 رولوراممد ۵۶ ۱۷۲60005 0110121 .1995 ۸0۵۸۰ -1 

۵ ۱/1۱ 60ا2طانامصا وم‌مهازه وومتع ۵۶ ومتاممت! ممتامنال۳۳۵ فقع فطلام۷۵0 .1993 ,اناومک .ل مه رب .ل مکلط0ظ -2 
71)04(:1041-۰ ,۹610060 امصصتصم ۵۶ امصیتول ,تب تفصتصصن 

۰ ,۳1۵2066 .ل 200 ولتطو2116 .2 مقصطه۱۱11 ۷۳۷ .ظ یه 920750 یک ر1(2۷165 یک .1 مهت)یازاا .ل همم .9 و8 ۸۷ مهمصهطانا -3 
0 0۲۵1116 «10)ما0م۳ج ققع تتقط) ۵۶ ممتامتنوعل 2 حطم و0عع تصقصتصصد که ممزاملهعم ۵۶ هن مطا مطتافصتاو۳۲ 
83:131-۰ مطمتالت ۱ ۵۶ لمصسمل صفتتبظ .عاممصه گم مموتنهم‌من :م۷ 1 

بو ت۷۵ وا م۷۷11 .عتدررامصه ممنووم‌جوه۲ 160امم۸ .1981 .اتمه ۲۱ صرح ریک .لا رکه1(]20 -4 

متصمعمصمطامه رها ممتامن0متم فقع ممتسممعصه عم فتاهتومموه مامصصلگ. خر ,1983 ,نا .۲ بو هه رب .ظ ۳600۲ -5 
.۰ -425 :(4)10 روع۲ا۱۵ «وماممطمع 1 لم)ممصمنمتص .فعلامها مهو صتوهت انم 

0 میامممعمه مه ومنامتصطمع) ۷۲۵0 ص1 0۶ ۳۲6۵۲ .2016 بعص2 ۷ 2۰ ۷۷۰ 20 رتهاولا۱6۸۵ ۸۵ .۲ ,220 با .۷ میک .2 ,۲۵ -6 
148-۰ :213 ,160001027 4ص ومممزمک ۳۵۵0 لقصصتصض .صمتاومون 1060 ج۵ وممر692 

تعح اجموعهاع0 اممنمد و ما ۵ ممتا‌ت0ع۳۳ .2008 بمصصت؟ ,۱۷ 4ص رتةزممتطه .9 مقلقمم96 ه و.ظ مصعصهاطن .-7 
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316-۰ :216 ,۱۷۱۵0611 ۳6۵۱۵02121 ۲۷607عوها۵ ,۷5 0۳0016160 ۵۶ 0۳60160160 

ص وام0مصه ممتامن0متن ققع ۷16۲0 م1 ۵ عو۲] ,1999 .12۵286 ,۲۱ .ظ 24 رلامعمطمو .ظ ولاه۲ .۲ .۸ و20955 ما .1 وبظ نک ما -16 
159)2(:145-1 ر5ع»ط10601۵ظ [2ع1)همجمط)۱۷۲2 .عمتامصتک امصتصنر 

کتاطاصهتمصه) . هدعم امه ۵۶ مماهه۷ ۵0 0۶ . ومتاولهامصعط. معقلله مهم نگ ._ ,2012 نک بل متصعطا2 -19 
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۵۵ وه آمنانه۲ ۵۶ مماقمصصتاز۴ .2011 .و2 .1 یه 24 رکلقلفصصفکا هر راقامطاصهن .0۵ رومهته۱ع۱۱ ۳۰ وی مصتطقه -20 
:۰ ,10۱21 16۲8181 زوعلنه[۳۵ تمصتهاه ۷ تممپزوهب نون وق .صمتاتل غمعصتصص؟ مر 0عون۲] دککنطاعلعع۳ محطم ۵۶ ومصرز] 
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027۷ ,مص1 مببااتاعما فمو ,8.0.1 مومعام ,۱1۵0005 10 وروی فقو م1 ,1999 ,9۸8 -21 

۶ لوصسمل ,حمتامنالمزج وفع مق 1۱ محم ممتاومعن معهای ۵۶ ومتامصتکا ,1994 ,۳۵۱۱ .۱۷ م۵ 220 الط .بط نگ وی م0001601٩‏ -22 


زابلی و همکاران,بررسی دقت برخی از مدل‌ها در برآورد ضرایب آزمون تولید گاز در سیلاژ ذرت ۳۷۷ 


۸۵1۴0۵1 5616066, 72: 29801۰ 

قوع 40صح متا باتااانافمون عصتفی تما فقو مهتم مصومو ۵۶ فمنلو پومنممه مطا که ممتاهمتصجماعل م1 .2002 ۳۲۰ رتمامه -23 
323-۰ :(153)5 ق۵۲۵0 ۷ ۱۷۲۵060186 0 میافک .اوعا 

فقو ۷۱۵۲0 ص11 عصلفدا فاصقصتصصن بق واتلهنان 1660 فصتصتصتهاهن1 ,2008 ۳۵۱1 ,۱ به مه ما10)مط90 .۲ و ما مت160650 -24 
عصتام۱۵0 جم ومطای/۷/۵ 4 1 نصا قفا ممت)ه)صمصصع؟ مامام‌مممه مه مصتعاتبظ .1 .منامتصطمعا ممتامنت2۳0۵ 
1 ,۴0۲2۵ 46 2ننال ر۵نامتصطهع) مملامنتلمج فقع فطا گم ممتههتاممرض :صمت)تابااظ کصفصتصیت۲ 

۶ ما0 .و0 10۲2۵ ۵۶ ممناوعع011 ۷۱/۲۵ م1 عط ۶0۴ مناوتصطهع) هع/۲۷0-9 ۸ .1963 .16۳۳7 ۸۵ با 4ص2 ریش ۷ .1 و1116 -25 
٩001607, 18: 104-7۰‏ 04صهاوجه طفتاترظ فطل 

مج تدای اصمعتمامل لحنامط وهای مامت 1۵۲ ولمطاع] ,1991 .12715 ۸ .ظ 224 رجمعاعحامط .ظ رل ول را۷20-۵065 -26 
3583-7۰ :74 رعممع۹01 نت۱ ۵۶ اقصیامژ .ممتا زانط تمصصتصح م6 م1612)10 م1 وم زمطم»۵۵1/۹2 ط‌تهاعط۵ظ 

لولعم ۵۶ ممتاه۷ 1۱۳۵۲ :ععتامصت! ممتام‌بالمم ققع ۷160 م1 عصتاع۷]۵0 ,2011 بصق1. مد .2 220 رعصع1 کر 9 و ,۷۷202 -27 
,137-0 :165 ,160000102 48ص ومصمزم ۳۵۵۵0 لمصصتصخ .8815مصظ ۵۶ ص0و1تهمرجمی هه فصمتاهناوع (ا) اوتامممممد۳۲ 

۶ صمته220ع0 ۵ ملهز امممتام مصا هن :2013 ,مممطمعظ ,را 20 رطع سنا 2 مها نک 9 مطللا .2 2 و ,۷۷202 -26 
920-۰ :7)6(۰ مأفطتمض .صمتاعممصصع؟ مل جا ترجه متقنونه م6 امعمصط (صنا) 21لاهمهممورم-15)10ع10 

5ص طمتاه)ممجصرع؟ امصتصصیا فا ومارععوع م6 واع0مصط 1هم)مصعطمصه ممو گم مصتاان1 ۵0۶ مموتنهم‌مصم) .2016 بطک1 رتل200 -29 
(صه1وع۳ م1) .35-47 :(5)3 رطه‌تهعوم؟ ممتامنل ۳۲۵ [محصتمد وقط ملظ عم میاموتصطهع؟ ممتامنا۳۳۵00 ققع ما مهصلل26001 

020۱۵121 6ظ) ۵۶ مصتام۱۷]۵0 ,1990 من ۷2۵ بک 24 روتتام‌طاصجمک ,۱۷ .۲ رتمم‌تنطامم‌مدمل بز وبط ۷ مع2۳10119 -30 
56)6(:1875-۰ راع۱۷۱۵۲۵01010 21)صممصمز صرظ صرح 0م1ا0جرخ .عنام 


نشریه پژوهشهای علوم دامی ایران ۵ م0۵۵۵ امصصنم۸ ۵۶ افصت ول صمتصرع رز 
جلد ۱۲. شماره ۰۴ زمستان ۱۳۹۹ ص. ۴۶۷-۴۷۹ رس 467-9 .۲ ,2021 ۱۷/۵۲ ,4 ,۷0 ,12 ,۷۵۱ 


00 6۵۵ 01 0۵0۵۵۵6۵ ۲۳۵۵ ۲۳۹)0۵۵۵)۵ و6 وز۱۷]۵0۵ ۵و ۵۶ مهت۸6 1۳۵ 
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"م۱۵۱6 جقفمم۱۷ مج کوصو اف مفی ر نا2200 ناف 
13-05-90 :1)60 وان 
21-04-00 :۸666060 


اقصتصی مصتاماويه »م1 0عونا. وهمطاعمه عط ۵۶ عم ود میمتصطهعا ممتام0متص فقع ۷۸۲۵ ۲7۱ صهمت6ام0۵ 1۱۸۳0۳۵0 
0صهوع۳ 15 طمت)2مابامصاً فطا هصتتیال 0عمن۵۳۵0 وهع ۵۶ محصتاام مطا رل‌مطاعصدر فنط ما .فتکنتو0عع گه متاممت! طمتاه/همصصررع] 
نومه ۵ ۵ )56 ماع مصتاا رها 0عصصرمهس ود مامت ممتام‌با0م0تم ققع ۵۶ ممتامتنوع0 آهم1امجهم‌طامه معط ,عبنه 2 و2 
,0۷76۷6۲ 0۳016 ممتام‌نا۳۵0 ققع عمصصتاوه 10 0676010۵60 معا ع2#ظ واع0مصظ تفممتممه لهتع۵7ز رزانمع»ع5 ,2۵061 
۶ ههام2 مطا مامهلاعع*ص1 ۲۵ و۷2 رنه فتطا ۵۶ حصتح معط رععملهت۵ظ1. ,طعاممه ماهتمعع )مط متع واع00ص عوعط) ۵۶ عمرمو 
0۰ 10۳۵ 2 ۶ »1امم! ممتاه)ممصع؟ اممتصند مصتاه01عتج عم امه عقعمتممه عمرمو 

۵ (11192فمه له درهل 60 همع 30 و0 مه فمامصصقع) و۹112 عهه رامه‌جصتنهمجه فطل ما 05مظ)عصه مصه واهزنه)۱۸ 
مج نا رصتهام0۳ عبت رهاظ متصحعنه) ممتندممرجم لههتصصمطه صرح تعنامم نا .عممتاوماناه ممتهاصهمصریع؟ قه ۷960 
۵ و1۳0 0مصتهاهان ف2/ ۱10 صمجصنر م1 ,عل‌مطامجه 0تهلصهاه مصلعن 0عصتصهاع0 ۳/۵۲۵ ومامصعه عطا ۵۶ (زناض 
کاععل؟ 2 مفصا 0عتعکفصه مج ۵۵۵۱60 ۱۷۵۲۵ 6تب؟ تفمتصت 0ما۵ع۱آهم م1 .8۵0 فهصتصمصه مطا مه فصصنهع 151212060 
د0 از تادنامنمتاومی وبا رطامامعومعطم فمه‌بوه1 عنام طوبامطا 0عتع)۱ وج متب معمصیاد م1 .بجماهتمطاق1 فطل مه 
معط که ممتانامی آهتممنمه 0ععگنه ۱۸16 0مرنمهر معط که 0ب ممجمییهر فظ1. بصمتهمانامصد مه 3960 21 2210161060 220 
مهم درم اصمعن0 معط مق فصن میهتهووو معط ما 0مامآمرمن مه میامتصطهع) ممتامنتلمعن وه ۷/۷(۰) 1:2 0۶ ۲۵010 
(ماماعطایه مه اما ولهز مب مه مهه‌تامتت صا 0متهمانامصا هتم ومامرصهه فص یت طمقع م1 .(قیرهل 6 60و12 صنط 
,72 ,48 ,24 ,20 ,16 م12 ,8 ,6 4 ,2 ,0 26 ۵۵50۲۵0 ۷۷28 0۳۵00۵۵0 225 0۶ مصصنیاآن۷ 1۳ .قاصهاها فطا قه 6۵88100760 ۱۷۵۲۵ 
۵ 210ها م1 ۵۶ ومحصت امن اج مصصیتام۷ فعع فطل ۵۶ ممتامت0هتهر م1 مممتاه‌مانمصد ۵ فتبامط 144 0ص 96,120 
0 1۳۵ ,100۷۵۵000 1۳۸ ۵۶ عجه 200 هه لصظ رد 0ووفمم0 ۱۸۵۲۵ کتالناوع۲ مرج وامعممط تهمصالممه عنام بوها 0عتهم‌جرمع 
7 200 (1۵0) ملاونع۱۵ (۳2۳۱) طمناطم۳۱ ور 2ظ۴) آهلامممممند۳۴ 060امص1 رقامصعمه لماممصتتممجع) وام2۵0 
۶ اوه ر(ظ۷۲۵) ۳0۲ ممهنانو صقفجد مصلوا 0ماقنلونه ۷۵۲۵ وام0مصهر فص ۵ 6 ۵۶ فععصهمع م1" ,(60۳۷) 
مطمتان00ع ص ۷۳۴(۰) تمه معهتهم‌هج معمصه 4صع (نتض/۵(٩)‏ عمتامتبعل متستاموواع صععص اعب‌نومر رح صمتاعمتصهامل 
مط ۶ ومبااه۷ 0۲۵010160 20 0۷۵0و۵0 ۵007۵۵۵) 227516 مم1ووم‌نوع1 تهمصتا مه ومع نار و( 12۱۷) ]وع) ۷۷/2090 -ظ1صانن رز 
عطا .ماج عطا عصتا ما 0615مصص عط ۵۶ برمهتتاهعه عطا ووعووه 10 ۱860 ۱۷۵۲۵ (فعصصلا ممتاهمانامصا کصمتعن که محصنتام۷ دهع 
۶ ووعص0همع هه (۲266 ممتامن0متم ققع مه مصصیتا۷۵ ققع 10ا0)حمصروه فطا) فتعامجصهتهم ممتاهتممع؟ تممتصب 0مامدتاوه 
«اماماموجصمی مصئعنا 2۵0وافصه م۱۷۵۲ (وعتاوتاهای ۱۷۲۳۳۴ قمع بام۳۳۸ خظ ر۵) عم طمعه حصمتق 0مصنداماه وعتعصصعتهم ]11 
101 1۲210101301260 

مج 30 ,0 چم (ه) محصیتام۷ وقع متامامصصره مصمتامن0۳60 ۵1 قصها ما مرمع مط 20 5لع0مصهر 0ع1بتای فط1 و6لاوع1 
۴ 98.37 0۶۲ هه معط ما ۷۵۲۵ 0615مصظ فص توها ۲0010060ض خر تعامجصحدم گم یله عطا لمح مصتالعمه ماه ورد 60 
40 1۱00018 1.00 صرح ۲2۳ م1۳ ,121۷۲ هد 200 بجر لمظر (60 بو جم ۵061 رن 6۵۲) 76.09 6۵ (0 وه چم 061مصه 002۷ 
۱۱۱ 0هتومصهن رممتناهعع 10 تفص مصتاهمت4صد رتم۲۵5 رقمباله۷ 12۲ همه ۱۹۴ اوعبسم۱ جع )ععطعنط عطا 
فص (اه۵007وع۲ ,2.879 220 2.591) 20015 ۳۸2۳۲۲ 24 000۷۲ 1 )۱۵۳/۵5 ۲28 مبتاج۷ ما۳۳ م1 .عام00ظ ۳۲۸2۲ 20 
۶ 11 0۲ ووه‌ص0ممع مطا عصتامناله7ه 1۵0۲ 1301602007 صرح وج 1560 15 میاه* باه۳۳۸ م1 ,(3.807) 200 ظعرظ ط1 )عولط و2 
0 4212 عصتا رز [0مصظ مطا ۵۶ تاج تمتتاهها ۵ فامموم۱0۵16 و(2870 ما تقو6۱۵) نه۴۵۷۸ ۵۶ ومیتاه۷ و1097 فط ,2۵016 
ص1 مبداه۷ ۱۷۲۳۴ م1 .0202 مصاان ما [06مددر ۲ صقطا پزممتنهعع تعطونط ه 20ظ واع0مجط ۲۸2۳۱ مه 00۷۲ روالناوع۲ وعطا ۵ 
0 136 ,۷۷۵۲۵5 قطان ص .(0.05ک) واع0مجط تعطان فط صقطا تعطونط زاامهم‌گتموزه ۳/۵۵ (5.527) 20061ظ ظع2ظ عطا 
معط صقطا 1۱0/۵۲ ۷۵۲۵ [06همجهر ظعرظ فطا ها (صمتد‌طانمصا ۵۶ فمصصتا اممعکن اج 0عمنمه وقع ۵۶ متام فط) فعباله۷ 
فص 1 ومتاعتاهاه ۱۱۱۷۷ 6 قااناوی: )۲25 (1۳۱۷۷) صمفع ۷۷ صزصاتتانا ی 0عفقظ .(0ممامعمصا قه۳ 1 ومیاله۷ 7۵0ر۵وها۵ 
0 ۲۸۳ )2 مصتاههم10 وراهبتامع۵و۲۵ ,0.675 200 0.705 ,0.691 ,0.392 ۱۷۵۲۵ ۴۵0015 1۷00 220 060۷ با۳2۳ ,۳۲۸۳ 
۶ 6)10ص1 ممناحاصهعصرع؟ معمصناد فطا متام ز0ع۳ن صا رراه۷ ۲۵5۵60۱ رزمدتتاهعع افعطونط مره ]1۵/۵5 مطا 20ظ واع0مصه 002 
۲۴ ۸ 10186 (< 3) صنا 107۷76۵1 عط) 20ظ [06مجدر زرط طا وعتیه فط 211 باقع صناد فا مه عصتت0عمظ مهقلته طتهع 


صق رصهعصه1! ررازوه/نصنا ممنه تلمبظ رمتانمته۸ ۵۶ ولنمد۲ راصهصانووه نز عمصعژمک لمحصتصظ ,0۳016950۲ اصهاک1وو۸ -1 
.2 مرحم حصق۱ ررازفه دنا همه تلصباظ متانه‌ته۸ ۵۶ واانموط رتمعمص‌تووعن1 مممممو: تمصصتصنظ رلعامنال ع2۳ ,1۷۲,۹6 -2 

( 7200110600860.260.11 :۴۱۵11 رمطاخظ مصلل ممموعتمن) -۳) 

101: 10.22067/1[251.71214 6 


زابلی و همکاران,بررسی دقت برخی از مدل‌ها در برآورد ضرایب آزمون تولید گاز در سیلاژ ذرت ۳۴۷۹ 


40 ۷6۲515 0090۳۷۵0 عطا معه ۵۵0 مملوومنوع۲ تههصل!ً م1 .کتاناوع؟ مطا هصاه‌ت0عتم صا [06ممدر کر ۵۶ همه 
1 م1 1 ۷10 5۱0۵ مه :0 اد )000تعص1 معع۵۲ه ممممعتنل اصمم‌نصوزه ۵ 60«مطو (قهامصصهتهم ط10ووعوع۲) وعبال۷2 
تحعصتا فطل مق 0مصنهاهاه مهتم ۶ ۵۶ وععصلممع فطا مق موه ۳۱۵۱۲۵۷7۵ ,(0.05>ع) ق0۵1مص 0160و فطل 
0۲00۱0۲10 وقع عط) ۵۶ ممتام01ععم تعااهها 2 20ظ 15ع0مصص 002۷ مه ۳2۲۱ ,1۲6۵۲695100 

مهفلنه جمه ۶ه متاممت! ممتامتمعصع؟ رمصصیر فط مصتامن0عتم طا پرمتناهعج 1۵9/۵۲ 20ظ 1ع0مصظ ظعرظ عط1 صمتمسلصمن 
0 20010108 صا 60وی و9 قاع0مصدر تمعصتاجمه قطان عقط 060همجصمرمعع1 و1 1۶ .0615مص 4160بتای مظان فظه ط 0عتهمعهع 
6 0۲۶ 0۶ ومتامص! ممتاهاصم‌جریع] [هصتصصیا فطل هصت)اهوتافع7ص1 عم عمج طع2ظ مطا 


۰ 2)10) مرج امصتصیت؟ رقاممصه تهممتامماظ بل‌مطاممط م۲2 2 ما ۵۶ 00000655 :۷۷۵۲۵8 مک 


